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Abstract.  The  Langmuir  plasma  probe  model  is  an  important  tool  in  spacecraft 
current  collection  and  charging  calculations.  In  ideal  geometries,  such  as  a  sphere  or  an 
infinitely  long  cylinder,  the  model  is  well  understood.  However,  the  realistic  geometries 
of  current  collectors,  or  spacecraft,  are  nonideal.  An  empirical  formula  for  a  Langmuir 
probe  with  a  given  nonideal  geometry  would  be  useful.  We  derive  such  a  formula  for 
the  SCATHA  satellite  by  using  the  SCIO  potential  data  obtained  during  electron  beam 
emissions.  The  satellite  rotated  perpendicular  to  sunlight  with  the  SC  10  booms  in  the 
equatorial  plane.  We  choose  one  special  mode  of  operation  during  a  quiet  space 
environment.  In  this  mode  the  beam  current  increased  continuously,  while  the  energy 
remained  constant.  We  analyzed  the  variations  of  the  vehicle  potential  responding  to 
the  unique  driving  factor,  the  beam  current.  To  provide  physical  explanation  to  the 
behavior  of  the  SC  10  potential  data,  we  model  the  interactions  between  the  beam, 
photoelectron,  and  ambient  currents.  We  present  an  algorithm  which  successfully 
yields  an  empirical  Langmuir  probe  formula  for  SCATHA,  from  which  we  obtain 
improved  estimate  of  ambient  electron  temperatures  and  densities.  The  results 
predicted  by  the  improved  Langmuir  probe  model  compare  favorably  with  the  very  few 
published  measurements  from  the  region. 


14.  SUBJECT  TERMS 

Spacecraft  charging.  Spacecraft 
SCATHA,  Satellite,  Geostationary 
beam 

Interactions,  Langmuir  probe, 

,  Geosynchronous,  Electron 

17.  SECURITY  CLASSIFICATION 

1  18.  SECURITY  CLASSIFICATION  I 

19.  SECURITY  CLASSIFICATION 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

NSN  7S4O-O1-28O-5S0O 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev  2-89) 

Pr«$cribM  by  ANSI  Std  139-^8 
29S-1Q2 


'*^^6  oixa 


Best 

Available 

Copy 


An  improved  Langmuir  probe  formuia  for  modeling  satellite 
interactions  with  near-geostationary  environment 

Shu  T.  Lai 


15  0  29 


JOURNAL  OF  GEOPHYSICAL  RESEARCH.  VOL.  99.  NO.  Al.  PAGES  459-467,  JANUARY  I,  1994 


An  improved  Langmuir  probe  formula  for  modeling  satellite 
interactions  with  near-geostationary  environment 

Shu  T.  Lai 

Phillips  Laboratory,  Hanscom  Air  Force  Base,  Massachusetts 


Abstract.  The  Langmuir  plasma  probe  model  is  an  important  tool  in  spacecraft 
current  collection  and  cha^'ng  calculations.  In  ideal  geometries,  such  as  a  sphere  or  an 
infinitely  long  cylinder,  the  model 's  well  understood.  However,  the  realistic  geometries 
of  current  collectors,  or  spacecraft,  are  nonideal.  An  empirical  formula  for  a  Langmuir 
probe  with  a  given  nonideal  geometry  would  be  useful.  We  derive  such  a  formula  for 
the  SCATHA  satellite  by  using  the  SCIO  potential  data  obtained  during  electron  beam 
emissions.  The  satellite  rotated  perpendicular  to  sunlight  with  the  SC  10  booms  in  the 
equatorial  plane.  We  choose  one  special  mode  of  operation  during  a  quiet  space 
environment.  In  this  mode  the  beam  current  increased  continuously,  while  the  energy 
remained  constant.  We  analyzed  the  variations  of  the  vehicle  potential  responding  to 
the  unique  driving  factor,  the  beam  current.  To  provide  physical  explanation  to  the 
behavior  of  the  SCIO  potential  data,  we  model  the  interactions  between  the  beam, 
photoelectron,  and  ambient  currents.  We  present  an  algorithm  which  successfully 
yields  an  empirical  Langmuir  probe  formula  for  SCATHA,  from  which  we  obtain 
improved  estimate  of  ambient  electron  temperatures  and  densities.  The  results 
predicted  by  the  improved  Langmuir  probe  model  compare  favorably  with  the  very  few 
published  measurements  from  the  region. 


1.  Introduction 

Ambient  electrons  and  ions  in  space  plasmas  impact  on 
the  spacecraft.  The  surface  potentials  come  to  equilibria  on 
timescales  of  milliseconds.  The  charging  time  is  so  short 
because  spacecraft  surface  capacitances  with  respect  to  the 
space  plasma  are  typically  small.  Some  dielectric  surfaces 
couple  to  subsurface  materials.  They  may  have  higher  ca¬ 
pacitances  and  accordingly  may  take  longer  amounts  of  time 
to  come  to  equilibria. 

At  equilibrium,  a  spacecraft  behaves  like  a  Langmuir 
plasma  probe  [Moit-Smith  and  Langmuir,  1926].  Since 
spacecraft  are  current  collectors,  their  potentials  are  gov¬ 
erned  by  the  balance  of  currents  as  given  in  the  Langmuir 
probe  equation.  The  control  of  a  spacecraft's  potential 
differs  from  that  of  a  laboratory  Langmuir  probe.  In  a 
laboratory,  currents  collected  by  Langmuir  probes  vary  in 
response  to  the  {q>plied  potential.  In  space,  however,  vehicle 
potentials  float  with  respect  to  the  ambient  plasma.  During 
charge-particle  beam  operations,  a  spacecraft's  potential 
varies  in  response  to  the  emitted  current  [Lai,  1989]. 

When  the  plasma  density  is  high,  the  current  density  is 
limited  by  the  Child-Langmuir's  law.  When  the  plasma 
density  is  low,  such  as  at  SCATHA  altitudes  (5-7  Rc),  the 
current  density  from  ambient  plasma  impacting  on  a  space¬ 
craft  surface  is  not  in  the  current-limiting  regime.  Instead, 
the  current  density  is  limited  by  the  orbital  angular  momenta 
of  the  incoming  charged  particles.  This,  we  refer  to  as  the 
orbit-limited  regime. 

In  the  orbit-limited  regime,  the  Langmuir  probe  equation 
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[MothSmith  and  Langmuir,  1926]  for  a  positively  charged 
probe  is 

mL(0)(i  -  «P  (-  ^)  =  0 

for  ed>  >  kT.  In  (I),  I ,{<!>)  and  /f(<^)  are  the  electron  and  ion 
currents,  respectively,  which  are  collected  from  the  local 
plasma  by  a  probe  charged  to  a  potential  relative  to  the 
plasma,  and  7  is  the  plasma  temperature.  In  (1),  -  -e, 

and  qi  -  e,  where  e  is  the  magnitude  of  the  elementary 
charge.  The  multiplicative  factor  p.  is  unity  if  the  probe  is  a 
perfect  sphere  and  2/7r'^  (-‘LI)  if  it  is  an  infinite  cylinder; 
the  power  factor  a  in  (1)  is  unity  and  1/2  respectively.  Since 
the  geometry  of  the  probe  is  neither  a  sphere  nor  an  infinitely 
long  cylinder,  n  and  a  have  neither  value. 

Lqframboise  arui  Parker  [1973]  generalize  the  concept  of 
orbit-limited  motion  to  include  spheroids.  They  conclude 
that  prolate  and  oblate  spheroids  have  an  orbit  limitation  in 
the  Laplace  limit  as  long  as  the  mqjor-to-minor  axis  ratios 
are  less  than  1.653  and  2.537,  respectively.  They  do  not  give 
a  values  for  specific  ratios.  Moreover,  SCATHA  is  not  a 
spheroid.  Thus  their  work  is  not  directly  applicable  here. 

Equation  (1)  describes  a  balance  in  which  the  total  current 
to  the  probe  equals  zero.  When  other  currents  are  involved, 
(1)  must  be  modified  accordingly.  For  example,  if  there  are 
currents  due  to  secondary  electrons  I,,  artificial  beam 
emissions  /«,  and  photoemissions  I  pi,  due  to  sunlight,  these 
terms  have  to  be  included: 

mL(0)(i  +  -I,  exp  =  /s  +  V  +  (2) 

The  magnetic  field  (~100  nT)  at  SCATHA  altitudes  is 
weak.  V  X  B  electric  fields  are  negligible  since  they  are  much 
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weaker  than  sheath  fields.  Also,  at  such  altitudes,  the 
ambient  ion  currents  are  typically  an  order  of  magnitude 
smaller  than  electron  currents  [Reagan  et  al.,  1981]  and  are 
unimportant  unless  the  spacecraft  is  charged  negatively  to 
thousands  of  volts. 

The  Langmuir  probe  equation  is  often  used  in  spacecraft- 
charging  calculations.  In  the  literature,  a  is  commonly  taken 
as  1  for  satellites  bearing  some  resemblance  to  spheres.  The 
main  objective  in  this  paper  is  to  obtain  an  improved 
Langmuir  probe  formula  for  the  SCATHA  satellite  by  deter¬ 
mining  the  value  of  a.  Toward  this  end,  we  must  first  identify 
and  estimate  the  various  currents. 

We  choose  to  study  the  SCATHA  satellite  measurements 
taken  in  the  period  2033  to  20S0  UT  on  March  II,  1981.  The 
altitude  was  42,600  km,  and  the  local  time  was  near  0500. 
This  day  is  chosen  because  of  the  following  reasons.  The 
ambient  environment  was  quiet  (not  stormy)  and  steady;  the 
24-hour  sum  of  3-hour  Kp  indices  was  14-  [Coffee,  1981]. 
During  the  period  considered,  an  electron  beam  was  emitted 
from  the  satellite.  The  beam  current  exceeded  the  ambient 
current  intercepted  by  the  satellite.  Therefore  the  spacecraft 
potential  was  controlled  by  the  electron  beam  emission. 
Furthermore,  the  beam  current  was  continuous,  and  the 
level  of  spacecraft  charging  was  positive,  in  the  range  0  to 
+ 170  V.  The  ion  current  can  be  ignored  under  these  circum¬ 
stances. 

In  section  2,  we  describe  the  SCATHA  satellite  and  the 
SCIO  instrumentation  on  board.  The  general  characteristics 
of  the  oscillations  of  SCIO  voltage  measurements  and  the 
observations  on  March  11,  1981,  are  described  in  section  3. 
To  study  in  detail  the  oscillation  behavior  on  March  1 1 , 1981 , 
we  present  our  theory  in  subsections  4.1,  4.2,  and  4.3.  In 
subsection  4.1,  we  model  the  photoelectron  currents  origi¬ 
nating  from  the  boom  surfaces  and  flowing  toward  the 
sateUite  body.  By  using  the  measurements  obtained,  we 
determine  the  photoemissivity  of  copper-beryllium,  the  sur¬ 
face  material  on  the  outer  segment  of  the  booms.  In  section 
4.2  we  delineate  three  interaction  regimes  depending  on  the 
contribution  of  various  currents  to  the  satellite  potential.  In 
section  4.3  we  choose  the  appropriate  regime  in  which  a  is 
the  only  parameter  to  be  determined,  and  we  present  an 
algorithm  for  calculating  a.  With  a  determined,  we  have  at 
hand  an  improved  Langmuir  probe  formula  for  SCATHA. 
Applying  this  formula  to  the  current-voltage  measurements, 
we  obtain  the  local  plasma  temperature  and  density.  Finally, 
in  section  5,  we  summarize  the  main  findings  of  this  paper 
and  compare  our  plasma  results  with  published  ones  ob¬ 
tained  with  different  instrumentation  and  techniques  under 
comparable  environmental  conditions. 

2.  SCATHA  Satellite 

SCATHA  was  launched  in  January  1979,  to  investigate 
natural  and  artificial  processes  controlling  charging  at  high 
altitudes.  Descriptions  of  the  experiments  on  SCATHA  are 
given  by  Fennell  [1982].  The  satellite  is  about  1  m  long  and 
1.6  m  in  diameter.  It  rotates  about  once  per  minute  with  a 
spin  axis  perpendicular  to  the  Sun-Earth  line.  SCATHA  is 
equipped  with  two  50-m  booms  (SCIO)  that  are  electrically 
isolated  from  the  satellite  ground.  The  surface  material  of  the 
outer  20  m  of  each  boom  is  made  of  an  exposed  copper 
beryllium  (CuBe)  wire.  The  inner  segment  is  coated  with 
kapton,  an  insulating  material.  The  SCIO  potential  repre¬ 


sents  the  difference  between  the  potential  ^uBc  of  the  tip  of 
a  boom  and  that  of  the  satellite  ground  [Lai  et  al.,  1986]. 
That  is, 

<t>  -  <^CuBe  ~  4>S-  (3) 

An  electron  beam  (SC4-1)  could  be  emitted  from  SCA¬ 
THA  with  various  energies  and  currents.  During  quiet  days 
in  sunlight,  SCATHA  normally  charges  positively  to  a  few 
volts.  The  emission  of  an  electron  beam  tends  to  raise  the 
satellite  potential  to  a  degree  that  depends  on  ambient 
conditions  as  well  as  beam  energies  and  currents.  When  the 
sateUite  rotates  in  sunlight,  the  amount  of  solar  illumination 
on  boom  surfaces  varies  sinusoidally,  as  does  the  photoelec¬ 
tron  current  from  the  booms  to  the  spacecraft. 

When  SCATHA  is  in  sunlight,  with  or  without  beam 
emission,  SCIO  potentials  show  oscillations  at  twice  the 
sateUite  rotation  frequency  [Lai  et  al.,  1986, 1987].  When  the 
sateUite  enters  eclipse,  the  amplitude  of  osciUation  decreases 
gradually;  this  evidence  supports  the  contention  that  the 
osciUation  is  due  to  the  effects  of  photoelectrons. 

When  SCATHA  charged  positively  as  a  result  of  electron- 
beam  emissions,  photoelectrons  from  the  booms  were  at¬ 
tracted  back  toward  the  satelUte  body  [Lai  et  al.,  1987].  The 
booms  form  part  of  the  sateUite  body's  electrical  environ¬ 
ment.  In  this  case,  the  satellite  body  not  only  interacts  with 
its  ambient  plasma  environment  but  also  with  the  booms 
(Figure  1). 

Because  of  the  high  secondary-emission  coefficient  of 
CuBe,  the  outer  sections  of  the  SCIO  booms  should  not 
charge  to  high  negative  potentials,  except  in  unusually 
energetic  plasma  environments  [Lai,  1991].  On  quiet  days, 
the  potential  </>cuBe  o*’  SCATHA  typically  varies  within  ±5V 
[Lai  et  al.,  1986].  When  the  spacecraft  potential  is  high 
compared  with  0cuBe>  ^  (3)  represents  a  good  approximation 
of  the  spacecraft  potential  with  the  sign  reversed,  that  is, 
d>  =»  There  are  several  instruments  for  measuring 

spacecraft  potential  on  SCATHA  [Fennell,  1982].  They 
confirm  that  SCIO  often  provided  good  approximate  mea¬ 
surements  of  the  sateUite  potential.  We  assume  that  SCIO 
measured  =  -  <^j ,  where  was  of  several  tens  of  volts 
and  was  controlled  by  electron  beam  emissions. 

3.  Observations 

We  describe  in  this  section  the  general  characteristics  of 
SCIO  potential  oscillations  and  the  specific  behavior  of  the 
oscillations  on  March  11,  1981. 

3.1.  General  Characteristics 

In  sunlight,  as  the  electron  beam  current  increases  from 
low  values,  the  spacecraft  potential  increases.  Furthermore, 
not  only  the  maxima  |  max((l>)  |  of  the  SCIO  potential  in  an 
osciUation  period  but  also  the  amplitude  |  max((/>)  -  min(<^)  | 
of  osciUation  increases.  The  extrema  of  the  oscillation  cor¬ 
relate  weU  with  the  Sun  angle  of  the  booms  [Lai  et  al.,  1987]. 
Minima  occur  at  0  =  0°  and  0  =  180°  and  maxima  at  90°  and 
270°.  Another  instrument,  SC2,  also  measured  the  potential 
of  the  spacecraft  body.  While  stUl  in  operation,  the  oscilla¬ 
tion  frequency  and  phase  of  SC2  potential  were  identical 
with  those  of  SCIO  during  electron  beam  emissions.  This 
indicates  that  the  oscillations  are  due  to  the  variations  of 

In  the  rest  of  this  section,  we  discuss  three  points  on  the 
SCIO  potential  data.  The  first  one  concerns  whether  the 
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Figure  1.  Schematic  diagram  of  the  SCATHA  satellite  with  beam  emission. 


boom  (SC2  or  SCIO)  potential  oscillations  always  depend  on 
vehicle  potential  modulation.  The  answer  to  this  question  is 
"no”  under  two  extreme  conditions  of  vehicle  charging.  In 
the  first  extreme,  the  vehicle  potential  relative  to  the  space 
plasma  is  a  few  volts,  and  is  nearly  constant,  as  during  quiet 
days  in  sunlight  without  beam  emission.  The  booms  may 
also  be  at  a  few  volts  relative  to  the  space  plasma.  Under  this 
condition,  there  is  a  twice  per  rotation  modulation  of  the 
boom  potentials.  Such  a  modulation  has  been  identified  to  be 
due  to  photoemission  from  the  booms  [Lai  et  al.,  1986; 
Craven  et  al.,  1987].  In  the  other  extreme,  the  space  plasma 
is  so  energetic  and  unusual  that  both  the  vehicle  and  the 
booms  are  charged  to  high  potentials  (a  few  kilovolts).  This 
condition  is  rare  [see  Lai,  1991a,  b]  and  did  not  occur  on 
quiet  days. 

The  second  point  to  discuss  concerns  the  use  of  SCIO 
rather  than  other  measurements  (SC2,  SC5,  and  SC9)  of 
vehicle  potential  on  SCATHA.  All  measurements  were 
approximate.  Because  of  the  short  length  (3  m)  of  its  boom, 
the  SC2  data  were  not  accurate  for  measuring  vehicle 
potential  relative  to  the  space  plasma.  The  SC2  instrument 
failed  eariy  in  1979.  The  SCS  data  below  about  100  eV  were 
inaccurate.  The  SC9  data,  which  measured  vehicle  potential 
by  identifying  the  shift  of  the  ambient  electron  distribution, 
were  also  inaccurate  at  low  potentials  and,  besides,  were 
sampled  at  a  slow  rate  (once  per  16  s). 

The  third  point  concerns  the  accuracy  of  SCIO  data  in 
representing  the  vehicle  potential.  At  high-vehicle  poten¬ 
tials,  the  Coulomb  potential  sheath  may  extend  far  beyond 
the  SO-m  booms,  the  data  would  merely  reflect  the  potential 
difference  between  the  vehicle  and  a  point  inside  the  sheath. 
For  potentials  below  about  300  V,  the  error  due  to  the 
Coulomb  sheath  is  negligible. 

3.2.  SateBtte  Potential  Osdllatioiis 

Figure  2  presents  data  acquired  on  March  11,  1981.  The 
beam  energy  was  constant  at  300  eV.  The  beam  current 


increased  continuously  from  near  zero  to  about  90  ftA. 
There  are  several  30-s  periods  of  calibration  dropouts.  Data 
taken  during  such  periods  are  ignored  in  our  study.  Oscilla¬ 
tions  in  the  potential  of  SCIO  detected  on  March  11,  1981, 
correlate  with  boom-Sun  angle  8,  with  spacecraft  potential 
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Figure  2.  March  1 1 , 1981 ,  measurements  on  SCATHA.  (a) 
SCIO  tb  potential  in  volts,  (6)  electron  beam  current  /{,  in 
microamperes,  and  (c)  SCIO  boom  Sun  angle  8  in  degrees 
are  presented  as  functions  of  time.  The  electron  beam  energy 
is  3M  eV.  The  dropouts  at  regular  intervals  are  for  calibra¬ 
tion. 
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nwxiiwa  occurring  wlwn  the  booms  were  parallel  or  antipar- 
altel  to  sunlight  direction  and  minima  when  the  booms  were 
perpendicular  to  it.  Starting  from  zero  beam  current,  the 
potential  oscillation  amplitude  increased  monotonically  with 
beam  cunent  until  a  critical  level  of  about  60  fiA  was 
reached.  The  amplitude  decreased  slightly  with  further  in¬ 
creases  in  beam  cunent. 

4.  Discusskm 

To  provide  a  physical  interpretation  of  these  observations, 
we  model  the  photoelectron  current  flowing  towards  the 
satellite  body  from  the  SC  10  booms  and  delineate  various 
regimes  of  interactions.  We  then  provide  an  algorithm  for 
determining  the  exponent  a  in  the  Langmuir  probe  equation 
(1). 

4.1.  Photodectron  Current  Modeling 

The  photoelectron  current  /ph(4i.  d)  from  the  booms  is  a 
function  of  the  Sun  angle  0.  Depending  on  the  potential  <!>,  of 
the  satellite  body,  some  fraction/of  this  current  flows  to  the 
main  body  of  the  spacecraft.  The  satellite  potential 
depends  in  a  self-consistent  maiuier  on  the  photoelectron 
current  /p(,(^,  6)  received  from  the  booms.  In  the  low- 
density  plasmas  at  SCATHA  altitudes,  the  orbit-limited 
Langmuir  plasma  probe  model  applies  for  the  collection  of 
ambient  currents.  The  current-balance  equation  for  the 
satellite  body  is 

M/,(0)(i  +  +  fph  («„  0)  =  /<,(<^,)  - 

(4) 


Figure  3.  Total  photoelectron  current  90°)  (equa¬ 

tion  (6))  arriving  at  the  spacecraft  body  as  a  function  of 
spacecraft  potential  4>a-  photoemissivity  ^pi,  used  is  3.S 
nA  cm“^  (equation  (5)). 


in  (5).  The  calculated  value  of  is  then  compared  with 
Kellogg’s  result. 

We  next  consider  the  photoelectrons  generated  on  a  boom 
surface  element  at  distance  r  from  the  main  body  of  the 
spacecraft.  Some  fraction  /  of  the  photoelectron  current 
moves  toward  the  satellite  with  the  rest  escaping.  The 
partition  of  photoelectron  current  depends  on  the  energies  of 
the  photoelectrons  and  the  sheath  potential  <Mr)  at  r.  In  this 
model  the  fraction /is  given  by  [Lai  et  al.,  1987) 


where 

fJOm 

fph  0)  =  2<f  I  dr/[<^(r)l/ph  sin  0  (5) 

Jo 

and  d  is  the  boom  diameter  /,(0)  is  the  ambient  current 
collected  if  the  spacecraft  potential  is  zero  with  no 
photoelectron  or  beam  emissions.  is  the  emitted  electron 
beam  current.  If  the  beam’s  energy  is  high  and  its  current 
density  low,  all  beam  particles  escape.  However,  if  the 
beam’s  energy  is  low  and  its  current  density  is  high,  some 
beam  particles  return  and  the  return  current  becomes 
nonzero. 

For  a  spherical  body,  the  power  a  of  the  orbit-limited 
current  collection  term  in  (1)  equals  unity;  for  an  infinite 
cylinder,  a  equals  1/2.  However,  the  SCATHA  satellite  is 
neither  a  sphere  nor  a  long  cylinder.  Rather,  it  has  a  short 
cylindrical  shape  with  nearly  the  same  length  and  diameter, 
llius  the  power  a  for  SCATHA  has  neither  value,  and  it  may 
fall  between  them. 

To  model  the  photoelectron  current  from  the  booms 
to  the  satellite,  we  assume  a  photoelectron  energy  spectrum 
and  a  sateUite  sheath  potential  profile  ^r)  as  a  function  of 
distance  r  from  the  satellite  surface.  Both  laboratory  and 
space  experiments  have  shown  that  a  Maxwellian  distribu¬ 
tion  is  a  good  approximation  for  describing  the  photoelec¬ 
tron  energy  spectrum  [Hinteregger  et  al.,  1965;  Whipple, 
1981).  The  photoemissivityjpt,  of  the  copper-beryllium  brom 
surface  material  on  a  rotating  satellite  has  been  estimated  to 
be  between  2  x  10“®  and  4  x  10”®  A  cm”^  [Kellogg,  1980). 
In  this  paper,  we  regard  jpi,  as  a  parameter  to  be  determined 


re<Kr) 

dE  E  txp  (-ElkTjt,) 

/[«^(r)]  =  -  (6) 

I  dEE  exp  (-EtkT^) 

Jo 

where  the  satellite  sheath  potential  <f>(r)  is  modeled  using  the 
Debye  form  [Whipple  et  al.,  1974): 

R 

<l>{r)  =  <^(0)  — — -  exp  (-r/Ao)  (7) 

r  "T  i\ 

Here  R  is  the  radius  of  the  satellite  body,  and  is  the 
Debye  distance.  For  SCATHA  environments,  we  assume 
that  the  Debye  length  Aj>  of  the  ambient  plasma  is  about  45 
m  [Aggson  et  al.,  1983],  and  the  photoelectron  temperature 
Tpi,  is  about  2  eV  [Whipple,  1981;  Lai  et  al.,  1986).  Using  this 
model,  we  have  computed  the  photoelectron  current  /ph(<^, 
90°)  going  toward  the  satellite  body.  The  results  appear  in 
Figure  3. 

The  maxima  and  minima  of  the  SCIO  potential  difference 
measurements  of  Figure  2,  plotted  in  Figure  4,  display  the 
emitted  beam  currents  as  functions  of  the  satellite  potential. 
Data  with  photoemission  (min  |  I  with  0  =  90°  or  270°)  from 
the  booms  are  plotted  in  Figure  5a,  and  those  without 
photoemissions  (max  |  0|  with  0  0  or  180°)  are  plotted  in 

Figure  5b.  At  low-beam  currents,  each  set  of  current-voltage 
data  varies  smoothly.  The  measurement  trend  suddenly 
deviates  from  the  trend  line  at  a  critical  current.  Without 
photoemission,  the  critical  current  is  about  60  /lA  with  the 
spacecraft  potential  at  about  220  V.  With  photoemission,  the 
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SC10  (Volts) 

Figure  4.  The  maxima  and  minima  of  SC  10  potential  4> 
(from  Figure  2)  as  functions  of  emitted  beam  current. 


critical  current  is  about  70  /xA  with  the  spacecraft  at  a 
potential  of  about  140  V. 

The  total  photoelectron  current  /ph(<^,  90°)  is  computed  for 
a  given  value  of  photoemissivity  jph  using  (5)  and  (6).  This 
computed  total  photoeiectron  current  /ph(<^>  90°)  is  then 
added  to  the  current  in  the  maximum  |  d)j  (B  =  0  or  180°) 
curve,  which  represents  the  spacecraft  body  potential  when 
there  is  no  photoeiectron  current  from  the  booms  to  the 


spacecraft  body.  The  result  of  this  addition  should  represent 
the  current  collection  of  the  satellite  body  during  maximum 
sunlit  conditions  at  the  booms,  that  is,  during  |  0|  min  (B  = 
90°  or  270°).  We  find  that  the  obtained  curve  fits  best  with  the 
experimental  data  points  of  min  |  <l>  \  when  the  value  of  the 
photoemissivity  jp^  (5)  is  about  3.5  x  10"’  A  cm"^  (Figure 
6).  This  value  of  jp^  determined  for  the  CuBe  surfaces  on  the 
SC  10  booms  of  SCATHA  agrees  with  that  (between  2  x 
10“’  and  4  x  10“’  A  cm“^)  estimated  by  Kellogg  [1980]  for 
the  CuBe  surfaces  on  the  Helios  spacecraft. 

4.2.  Interaction  Regimes 

The  satellite  potential  (^,(0)  oscillates  as  the  satellite  and 
booms  rotate  in  sunlight.  The  amplitude  of  the  potential 
oscillation  is  given  by  Sif>  =  <^,(0  =  0  or  180°)  -  <^,(0  = 
90°  or  270°).  As  the  beam  current  /*  increases,  so  does  the 
When  the  satellite  charges  to  near  beam  potential,  the 
beam  current  does  not  escape  completely  [Olsen,  1989]. 

Three  regimes  of  interactions  1,2,  and  3  can  be  identified 
in  cases  with  and  without  photoemission  (Figures  5a  and 
5b).  In  regime  1  and  2,  the  /p  is  unsaturated,  and  the  return 
current  is  zero.  In  regime  1 ,  some  photoelectrons  from  the 
satellite  body  escape  because  the  satellite  positive  potential 
4>s  is  low.  The  amount  of  photoeiectron  current  leaving  the 
satellite  body  is  a  function  of  satellite  potential  </>s  ■  In  our 
model  the  current  balance  equation  is 

e<f>Y 

1+^1  +Ipi,(4>s>6)  =  I„(il>,)  +  Jpi,(4>„e) 

(8) 


I  I  m  I  I  m 


Figure  5.  (a)  The  emitted  beam  current  and  the  SCIO  potential  with  photoemission  (minimum  |  <]>\ 
with  0  =  90°  or  270°).  (b)  The  emitted  beam  current  /j,  and  the  SCIO  potential  without  photoemission 
(maximum  |  <l>  \  with  0  =  0  or  180°).  The  three  regimes  in  each  case.  Figure  5a  or  5b,  are  discussed  in  the 
text.  The  cartoons  below  the  x  axis  show  schematically  the  physical  processes  in  the  regimes. 
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Figure  6.  Fits  to  determine  the  photoelectron  current.  The 
lowest  curve  is  a  polynomial  fit.  The  difference  between  the 
two  data  sets  shown  in  Figure  4  is  fitted  by  using  (5).  The 
upper  solid  curve  is  obtained  by  using  =  45  m  and  = 
3.5  nA  cm~^  and  the  dash  curve  Aq  =  12  m  and  =  6.5nA 
cm~^. 


where  the  photoelectron  current  leaving  the 

satellite  body. 

In  regime  2,  the  photoelectron  current  from  the  satellite 
body  is  very  small  because  the  satellite  body  potential 
energy  e<l>,  is  high  compared  with  the  photoelectron  energy: 


filAO) 


+  /ph(0..  9)=^tA<t>s) 


(9) 


By  extrapolation  [Gonfalone  el  al.,  1979],  the  regime  2  curve 
intercepts  the  y  axis  =  0)  at  about  10  fiA  (Figure  5b). 
This  determines  fi/,(0)  to  be  approximately  10  itA  in  (4) 
which  gives  /j,(0)  =  /a/,(0)  =  10  /iA  at  the  intercept.  Using 
the  surface  area  D  =  9.05  m^  of  SCATHA  and  the  average 
ambient  current  density  =  0.115  ±  0.10  nA  cm"^ 
obtained  in  a  45  days  average  in  1979  at  SCATHA  altitudes 
[Purvis  et  al.,  1984],  one  obtains  a  result  for  the  average 
current  intercepted  by  SCATHA  as  07,  =  10.4  ±  9.1  x 
10"*  A  if  At  =  1  or  9.5  ±  8.3  X  10"*  A  if  p,  =  1.1.  This 
average  result  is  of  the  same  order  of  magnitude  as  the  7,(0) 
determined  («10  /uA),  the  ambient  current  intercepted  by 
SCATHA  on  March  11,  1981. 

In  regime  1,  the  data  points  deviate  from  the  curve 
extrapolated  from  regime  2,  because  of  photoelectrons  leav¬ 
ing  the  spacecraft  body.  On  the  other  hand,  very  few  of  the 
photoelectrons  from  the  booms  can  reach  the  spacecraft 
body,  because  of  the  low  attraction  offered  by  .  In  regime 
2,  a  multibody  interaction  occurs  between  the  satellite  body, 
the  booms,  the  electron  beam,  and  the  ambient  plasma.  This 
regime  will  be  discussed  further  in  section  4.3.  In  regime  3, 
beam  saturation  occurs,  and  part  of  the  beam  current  re¬ 
turns;  this  regime  is  not  of  main  interest  here.  Physical 
processes  characteristic  of  the  three  regimes  are  shown 
schematically  at  the  bottoms  of  Figures  5a  and  5b. 


(VOLTS) 

Figure  7.  Solutions  to  (12;  to  determine  the  parameter  a. 


4.3.  An  Algorithm  for  Determining  a 
We  now  offer  an  algorithm  for  determining  the  value  of 
power  factor  a  (1).  When  there  is  no  photoemission  in 
regime  2  (Figure  5i>),  the  only  currents  involved  are  the 
ambient  electron  and  ion  currents  and  the  beam  current 
Ih .  Secondary  electrons  have  insufficient  energy  to  escape, 
since  the  <^,  >  +30  V.  Furthermore,  neither  photoelec¬ 
tron  current  /ph  nor  beam  return  current  I,  is  involved 
in  this  regime.  The  beam  current  /*  provides  the  driving 
force;  the  response  <l>  is  then  a  function  of  /{,  only.  Unlike 
/ph  and  I,,  the  beam  current  is  controlled  and  there¬ 
fore  known. 


(10) 


Differentiating  (10)  by  e4>  gives 


afi  ( 

/'.(.0)~MO)  1+- 


where  /J,  denotes  differentiation  of  7*  by  e<f>.  Dividing  (10) 
by  (1 1),  we  obtain 


lb  1 

-  =  -(Ar+c<fr).  (12) 

lb  “ 

Equation  (12)  is  a  simple  formula  enabling  a  and  AT  to  be 
determined.  If  the  data  7h(<^)  satisfy  the  Langmuir  probe 
equation  (10),  they  should  satisfy  (12).  An  algorithm  based 
on  (12)  is  as  follows.  Plotting  the  data  7^/7),  as  a  function 
(ordinate)  of  e(t>  (abscissa)  should  yield  a  straight  line  with 
I/a  as  the  slope.  Once  a  is  determined,  the  ambient  electron 
temperature  T  can  be  determined  from  the  intercept  (at  = 
0)  which  equals  kHa. 

Using  data  from  Figure  5b,  Figure  7  shows  a  plot  of  IbH'b 
as  a  function  of  e<f>.  In  Figure  7  the  data  7^/7),  falls  into  a 
straight  line  for  the  potential  range  from  ift  «=  25-160  V.  This 
linear  dependence  of  7^/7],  on  e<i>  is  as  predicted  by  (12). 
Outside  this  range,  the  current- voltage  behavior  deviates 
from  the  simple  Langmuir  probe  equation  given  by  (10). 
Below  0  20  V,  some  secondary  electrons  and  photoelec- 
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Irons  from  the  spacecraft  body  escape,  and  above  =  160 
V,  some  beam  electrons  return.  Between  the  two  critical 
values  of  potential  0,  //,//'»  falls  on  a  straight  line,  as 
predicted  by  (12).  Prom  the  slope  of  this  line,  we  obtain  the 
power  factor  a  =  0.774.  This  fitted  value  of  a  lies  about 
halfway  between  the  values  for  a  sphere  and  an  infinitely 
long  cylinder. 

After  obtaining  a,  one  can  obtain  the  ambient  electron 
temperature  AT  from  the  intercept  in  Figure  7.  The  intercept 
equals  AT/a  and  gives  AT  =  23.2  eV.  This  value  is  compa¬ 
rable  to  AT  =»  64  eV  [Whipple,  1981]  which  was  measured  by 
a  different  method  on  the  ATS-S  satellite.  Since  the  environ¬ 
ment  (e.g.,  plasma  sheet  and  plasmasphere)  is  unknown,  the 
comparison  is  not  meant  to  be  strict.  We  merely  note  that  the 
measurements  are  of  the  same  order  of  magnitude. 

Taking  AT  ■'=  23.2  eV  and  the  ambient  current  /,(<^  =  0) 
»  10  fiA  obtained  from  the  intercept  (section  5),  we  can 
determine  the  ambient  electron  density  n,  by  the  following 
equation: 

,  en,D  (SAT\ 

/,(0)  =  i  =  — —  -  «10mA  (13) 

*  4  \irm  j 

which  gives  n,  =  8.47  cm“^  for  /i  =  I  and  7.70  for  =  1. 1. 
We  cannot  determine  p  without  knowing  7,(0)  indepen¬ 
dently. 


5.  Summary  and  Conclusions 

The  equilibrium  potential  of  a  spacecraft  is  governed  by 
current  balance  as  given  in  the  Langmuir  probe  equation. 
The  power  a  in  the  Langmuir  probe  equation  [Mott-Smith 
and  Langmuir,  1926]  is  well  known  to  be  I  or  1/2  for  a 
spherical  or  infinitely  long  cylindrical  probe,  respectively. 
For  most  spacecraft,  the  geometry  is  neither  a  sphere  nor  an 
infinite  cylinder.  For  improved  spacecraft  charging  model¬ 
ing,  it  is  better  to  determine  and  use  the  correct  value  of  a 
rather  than  taking  a  =  I  as  commonly  practiced.  This  leads 
to  a  lower  estimate  of  AT  and  a  higher  estimate  of  the  plasma 
density. 

With  beam  emissions,  the  beam  current  emitted  is  likely 
known  except  when  there  is  beam  return.  From  current- 
potential  measurements,  it  is  then  possible  to  deduce  the 
value  of  a.  Before  one  can  apply  an  algorithm  to  deduce  the 
power  a,  one  has  to  identify  and  disentangle  all  the  various 
interactions  between  the  spacecraft  body,  the  booms,  the 
ambient  plasma,  and  the  beam  emitted  or  even  returned.  We 
have  taken  the  March  1 1,  1981,  current-voltage  measure¬ 
ments  obtained  on  SCATHA  satellite  for  a  case  study.  The 
day  was  chosen  because  there  was  no  storm  and  the  period 
chosen  was  the  only  one  in  which  the  SCATHA  electron 
beam  current  was  increasing  continuously  instead  of  chang¬ 
ing  in  large  steps. 

5.1.  Key  Observations 

The  SC  10  potential  difference  measurements  oscillated  at 
twice  the  satellite  frequency  with  their  maxima  occurring 
when  the  booms  were  parallel  or  antiparallel  to  sunlight  (0  = 
0°  or  180°)  and  the  minima  occurring  when  the  booms  were 
perpendicular.  The  amplitude  of  oscillation  increased  with 
the  beam  current  until  the  satellite  potential  was  near  the 
beam  energy. 


5.2.  Key  Interpretations 

When  the  booms  were  aligned  with  sunlight,  there  was  no 
photoemission  emitted  from  the  boom  surfaces.  When  the 
booms  were  perpendicular  to  sunlight,  the  photoemission 
from  the  booms  was  maximum.  As  the  emitted  beam  current 
increased,  the  spacecraft  potential  increased.  As  a  result,  the 
spacecraft  sheath  engulfed  part  of  the  booms,  and  therefore 
some  photoelectrons  from  the  booms  were  attracted  toward 
the  satellite  body.  The  amount  of  photoelectron  current,  and 
therefore  the  amplitude  of  oscillation,  increased  with  the 
spacecraft  potential.  When  the  spacecraft  potential  was  near 
the  beam  energy,  the  beam  partially  returned.  We  have 
identified  three  different  interaction  regimes  (Figure  5).  In 
regime  1,  photoelectrons  were  emitted  from  the  spacecraft 
body.  Also,  in  this  regime,  the  approximation  <t)  =  -(t>s  nt^y 
not  be  accurate.  In  regime  3,  partial  beam  return  occurred. 
In  between  these  two  regimes,  we  have  modeled  the  inter¬ 
actions. 

5.3.  Calculations  and  Results 

1.  Using  a  simple  photoelectron  current  partition  model 
to  fit  the  measured  oscillation  amplitude,  we  found  that  the 
photoemissivity  7ph  of  the  boom  surface  was  about  3.5  x 
10"’  A  cm“^. 

2.  The  ambient  current  7(0)  can  be  deduced  by  extrapo¬ 
lating  the  regime  2  curve  to  intercept  the  =  0  axis.  From 
the  intercept,  we  obtained  the  ambient  current  7(0)  inter¬ 
cepted  by  SCATHA  as  10  pA. 

3.  For  Langmuir  probe  modeling,  we  chose  the  current- 
voltage  measurements  in  regime  2,  when  the  booms  were 
aligned  with  sunlight.  No  photoemission  from  the  booms 
was  involved  with  these  measurements.  No  partial  beam 
return  occurred  in  this  regime.  The  spacecraft  potential  <i> 
was  a  function  of  the  beam  current  I/, ,  the  ambient  current 
7(0),  and  the  ambient  temperature  T^  only  (10).  We  plotted 
7j,/7'fc  as  a  function  of  <l>.  The  slope  gives  a  and  the  intercept 
ATJa.  With  a,  AT^,  and  7(0)  known,  we  can  deduce  the 
ambient  electron  density  n,.  The  results  obtained  are  a  = 
0.774,  ATf  =  23.2  eV,  and  n^  =  8.47  cm"^  for  =  1  or/?, 
=  7.70  for /it=  1.1. 

5.4.  Comparisons  With  Other  Measurements 

The  calculated  photoemissivity  result  agrees  with  the 
measurements  by  Kellogg  [1990].  The  ambient  current  value 
obtained  agrees  with  the  statistical  measurements  of  Purvis 
et  at.  [1984].  It  is  interesting  that  the  value  of  a  =  0.774  lies 
between  the  two  known  values  1  and  1/2,  corresponding  to  a 
perfect  sphere  and  an  infinite  cylinder,  respectively.  Since 
the  geometry  of  SCATHA  resembles  a  short  cylinder,  the  a 
result  seems  reasonable.  The  ambient  temperature  IcT,  is 
comparable  to  the  value  AT^  =  64  eV  obtained  on  ATS-5 
[Whipple,  1981]  at  the  geosynchronous  environment  by 
mtans  of  different  instrumentation  and  technique.  Since  the 
environment  is  unknown,  we  merely  note  that  our  result, 
ATf  =  23  eV,  is  not  unreasonable  for  near  geosynchronous 
orbit,  on  a  quiet  day  near  local  dawn,  at  5.5  (perigee  for 
SCATHA). 

We  compare  our  /?,  result  with  the  statistical  results  of 
GEOS  measurements  obtained  by  means  of  relaxation 
sounding  at  6.6  7?^  [Higel  and  Lei,  1984].  During  the  p>eriod 
considered,  the  altitude  of  SCATHA  was  approximately 
40000  km  (=7.3  R^)  at  0500  LT  and  the  IKp  was  14- .  At 


m 
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the  values  <rfXlkp,  the  GEOS  plasoia  densities  at  6.6 
Rg  and  QSOQ  LT  were  about  3  to  7  cin~^  for  IKp  =  13.  1  to 
2  cm'^  for  SJCp  »  IS,  and  8  to  9  cm'^  for  IKp  =  21  [see 
Higel  and  Lei,  1984,  Figures  3E,  3A,  3B].  In  t^  compari¬ 
son.  our  local  plasma  density  n,  is  on  the  high  side. 

$.5.  Finai  Conraienis 


V  = 


—  z 
kT 


0 


(17) 


For  J  =  0,  p(0)  =  ir'^^/2.  For  large  s,  g(s)  -*  l/2s,  one 
approximates  in  (16)  by  l/2n'^.  then  approximates  1 

+  l/2i)  by  (1  +  Tj)*'^  and  obtains 


We  have  assumed  a  Debye  length  Aq  =  4S  m  [Aggson  ei 
al.,  1983]  to  fit  for  the  photoelectron  current  (top  solid  line  in 
Figure  6).  In  retrospect,  we  use  A|,  =  12  m  (kf,  »  23.2  eV 
and  n,  »  7  cm  ~^)  to  fit  (dashed  line  in  Figure  6),  yieldingjp), 
••  6.5  X  10"®  amp  cm~^.  This  photoemissivity  value  is 
higher  than  Kellogg's  [1980]  estimate. 

In  our  algmithm  (12)  we  have  avoided  any  photoelectron 
current  flowing  from  the  booms  to  the  vehicle  body  by 
selecting  the  branch  of  SC  10  potential  data  points  (Figure  4) 
with  Sun  angle  0  =  0°  or  180°.  Any  error  in  the  initial 
estimate  of  plasma  density,  or  Debye  length,  would  not 
affect  the  algorithm  for  determining  1,(0),  a,  kT,,  and  n,. 

Finally,  we  briefly  comment  that  the  presence  of  a  satellite 
may  affect  its  local  plasma  density.  Enhancement  of  local 
plasma  density  has  been  observed  by  Olsen  et  al.  [1981], 
who  attributed  the  cause  to  possible  potential  barriers  due  to 
differential  charging  on  the  satellite  surface.  In  our  case, 
electron  beam  emissions  may  also  be  plausible  reasons.  For 
purposes  of  designing  better  experiments  for  improved  de¬ 
termination  of  Langmuir  probe  parameters  for  spacecraft 
and  the  parameters  of  the  ambient  plasma,  we  suggest 
emitting  a  high-energy  (several  keV)  electron  beam  at  night. 
High-energy  beam  would  generate  less  electrons  in  the 
vicinity  of  the  spacecraft.  Without  photoemission  from  the 
spacecraft,  the  Langmuir  probe  current  at  =  0  reduces  to 
7,(0)  regardless  of  the  geometry,  sphere,  or  cylinder.  This 
would  determine  the  ambient  current  without  extrapolation 
or  the  ambiguity  factor  p.  between  1  and  1.1. 

Appendix 

This  appendix  discusses  the  differences  in  functional 
forms  of  Langmuir  probe  current  due  to  spherical  and 
cylindrical  geometries.  The  flux  J  collected  by  a  spherical 
probe  is  given  by  [Mott-Smith  and  Langmuir,  1926] 

where  is  the  random  flux  n.„(kTI2irm)''^ ,  n*  is  the 
plasma  density  at  infinity,  and  m  is  electron  mass.  The  flux 
J  is  related  to  the  current  /  by  a  multiplicative  factor  of 
surface  area. 

For  a  cylindrical  probe,  the  flux  J  collected  is  given  by 
[Mott-Smith  and  Langmuir,  1926;  Laframboise  and  Parker, 
1973] 

y  =  7„(2/7r''^)[Tj''^+p(i7''^)]  (15) 

where 

g{s)  =  J  exp  (s^)  erfc  (s) 

=  exp  (j^)  f  exp  (-/^)  dt  (16) 


J 


(18) 


The  functions  [ij’^^  +  P(ij'^^)l  and  (1  +  tj)*^  versus  tj 
[Mott-Smith  and  Langmuir,  1926;  Swift  and  Schwar,  1977] 
already  show  little  difference  for  >  kT  and  no  visible 
difference  at  all  for  e^  >  IkT.  This  property  further 
supports  that  our  starting  point  of  the  data  fitting  (section 
4.3)  at  about  30  eV  (see  Figure  7),  below  which  the  data 
deviate  from  the  Langmuir  probe  form  considered. 
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